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Abstract 

We study a scalar field theory with a non-standard kinetic term mini- 
mally coupled to gravity We establish the existence of compact boson stars, 
that is, static solutions with compact support of the full system with self- 
gravitation taken into account. Concretely, there exist two types of solutions, 
namely compact balls on the one hand, and compact shells on the other hand. 
The compact balls have a naked singularity at the center. The inner boundary 
of the compact shells is singular, as well, but it is, at the same time, a Killing 
horizon. These singular, compact shells therefore resemble black holes. 
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1 Introduction 



This paper investigates self-gravitating compact solutions of a non-linear scalar 
field theory with a non-canonical kinetic term. Recently it has been established 
that relativistic non-linear field theories may have static (solitonic) solutions with 
compact support, such that the fields take their vacuum values identically outside 
a compact region. Solitons of this type are called "compactons". At the moment, 
there are two known classes of field theories which may have compacton solu- 
tions. One possibility is that the (scalar) field of the theory has a potential with 
a non-continuous first derivative at (some of) its minima, a so-called V-shaped 
potential |jT| - [[TOll . The other possibility consists in a non-standard kinetic term 
(higher than second powers of the first derivatives) in the Lagrangian [fTTI . [[T2l . a 
so-called K field theory. It may be interesting to mention at this point that K field 
theories have found some applications already in cosmology as a candidate for 
dark energy, where theories of this type are known as "K essence", or "general- 
ized dynamics", see, e.g., ||T3l - [fTSll . They have been also applied to topological 
defect formation, see, e.g., |fT9ll , [|20ll . 

In both cases of V-shaped potentials and of K theories, respectively, the expo- 
nential approach to the vacuum typical for conventional solitons is replaced by a 
power-like approach. Further, the vacuum value is reached at a finite distance (the 
"compacton boundary"), and the second derivative of the field is non-continuous 
at the compacton boundary. In some cases, the compacton solutions are, therefore, 
weak solutions (i.e., they do not solve the field equations at the boundary) which, 
in these cases, is not problematic, because the space of weak solutions is the ade- 
quate solution space for the corresponding variational problem. When gravitation 
is included, the stronger continuity requirements for a space-time manifold might 
render such a behaviour problematic, so let us emphasize already at this point that 
this does not happen for the self-gravitating compact solutions considered in this 
paper, i.e., all solutions are strong solutions of the Einstein equations. 

Like in the case of conventional solitons, for compactons it is also much sim- 
pler to find solutions in theories in l-i-l dimensions. We remark that for some 
non-relativistic theories (generalizations of the KdV equations) compactons were 
already found in fTP\, fTT\. The most direct generalization of topological com- 
pactons to higher dimensions faces the same problems like in the case of con- 
ventional solitons, and also one possible solution is the same, namely the in- 
troduction of additional gauge fields, like in the case of vortices and monopoles 
||23l . Other possibilities to find higher-dimensional compactons consist in choos- 
ing more complicated topologies for the target space fi24il or in allowing for a 
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simple time dependence in the form of Q balls [|6||, Q. 

A further natural question to be asked is whether these compactons may be 
coupled to gravity and whether such self-gravitating compactons lead to some in- 
teresting applications. A first application has been to brane cosmology, where 
the self-gravitating compactons provide a model for a thick brane with a strictly 
finite extension, as well as an automatic confinement of all matter fields to the 
brane [|25l , [|26l , [|27l . As this thick brane is of the domain wall type, the resulting 
system is essentially a one-dimensional problem. The coupling of the 3-1-1 di- 
mensional compact Q balls of the V-shaped class of models of [TJ to gravity was 
recently studied in [|28ll mainly using numerical methods. 

It is the main purpose of the present paper to study in detail the coupling of 
a 3-1-1 dimensional radially symmetric compacton of a K field theory (with non- 
standard kinetic term) to gravity. For solutions to a scalar field plus gravity theory 
the notion of a "boson star" has become standard in the literature, therefore we 
will call generic solutions of our model "compact boson stars", and reserve more 
specific terms ("compact balls", or "singular shells") for the different types of 
solutions we shall encounter. For boson stars with a conventional kinetic term of 
the scalar field, there already exists a large amount of literature. Here one first 
criterion for a classification is whether the scalar field is real or complex. In the 
case of a real scalar field, a static, spherically symmetric metric requires a static, 
spherically symmetric scalar field. The resulting field equations are known not to 
support solutions when gravity is neglected, because their existence is forbidden 
by the Derrick theorem. For the corresponding system with gravity included, finite 
(asymptotic) mass solutions do exist, but they are beset by naked singularities, 
see e.g. [291 - [|32| . For a complex scalar field, the Q-ball type ansatz </)(x) = 
exp(— ico't)/(r) still is compatible with a static, spherically symmetric metric. In 
addition, this ansatz allows for regular, finite energy solutions already in the case 
without gravity. As a consequence, regular finite mass solutions for the full system 
with gravity included do exist and have been widely studied, see e.g. [33J - [|35ll . 
Reviews on boson stars may be found, e.g., in [l36l - [l38l . We remark that the 
system studied in the present paper is different in this respect, because, due to the 
non-standard kinetic term, static finite energy solutions in the non-gravitational 
case are not excluded by the Derrick theorem, and they do indeed exist, as we 
shall see in the next section. 

Our paper is organized as follows. In a first step, in Section 2 we introduce 
the simplest K field theory which gives rise to non-topological compact balls in 
flat Minkowski space (i.e., without gravity). We discuss general features of the 
resulting compact ball solutions and calculate some explicit solutions by numer- 



2 



ically integrating the ODE for spherically symmetric compact balls. In Section 
3, we couple the model of the previous section to gravity via the usual, minimal 
coupling. We derive the Einstein equations for spherically symmetric configura- 
tions and discuss in detail their properties. We then present some explicit compact 
boson star solutions by numerically integrating the ODEs for the spherically sym- 
metric configurations, and, finally, we investigate the singularities which form in 
these solutions. We find, in fact, two types of solutions which have two different 
types of singularities. The first type has naked, point-like singularities and is, thus, 
consistent with the "no scalar hair" conjecture. The second type has a surface-like 
singularity that is, at the same time, a singular boundary of space-time, and the lo- 
cus of a (Killing) horizon. This second type of solution, therefore, exhibits some 
similarity with a black hole. Section 4 contains a discussion of these results as 
well as some speculations about possible astrophysical or cosmological applica- 
tions of compact boson stars of the type studied in this paper. 



2 The non-gravitational case 

We want to study the simplest possible scalar field theory in 3-1- 1 dimensions with 
a quartic kinetic term and a standard potential which gives rise to the formation 
of non-topological compact solitons (i.e. compactons). Therefore we choose the 
Lagrangian 

C^-X\X\-V{<t>), (1) 

where 

X = \g^^d^(t>d,(t>, V{(t>) = \(t>\ (2) 

We remark that, in the case at hand, the Derrick scaling stability criterion does 
not exclude the existence of static solutions, thanks to the non-standard, quartic 
kinetic term. The absolute value symbol in the kinetic term is necessary to ensure 
boundedness from below of the energy. For static configurations this absolute 
value symbol is immaterial in the case without gravity, because the metric func- 
tions are given functions with a fixed sign. For the case with gravitation included 
this is not guaranteed a priori because the metric functions have to be determined 
from the Einstein equations and could, in principle, change sign (as happens, e.g., 
in the case of the Schwarzschild solution). We will find, however, that such a sign 
change is not possible in our model. 
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The line element ds^ of flat spacetime in spherical polar coordinates reads 

ds^ = -dt^ + cir^ + r2(d^2 ^ gi^2 Q^2,^ 

Now we assume spherical symmetry = 0(r) then the Euler-Lagrange equation 
takes the form 

(r _ 2A0 = 0. (4) 

It is convenient to introduce the new variable s = r^/^ instead of r. In the variable 
s equation © takes the form 

(f)'^(j)" - 54AsV = 0. (5) 

The advantage of the new variable s is that solutions have a power series expansion 
about zero in the variable s, but not in r. We remark that this will no longer hold 
in the case with gravitation included. There solutions will have a regular power 
series expansion about r = in the variable r. For brevity, in the sequel we shall 
call the variable s "radius", as well, and distinguish the two variables only by their 
respective letters Q 



2.1 Expansion at the boundary 

As always in the case of compactons, we now assume that there exists a certain 
radius s = S where the field takes its vacuum value, and its first derivative van- 
ishes. Further, the second derivative is nonzero from below (for s < S), whereas 
it is zero for s > S (such that takes its vacuum value zero for s > S). Indeed, 
expanding around s = S, 

(l>is)=Y.f,is-Sy for s<S (6) 

j 

and assuming that 0(s = S) = (f)'{s = S) = we get a cubic equation for /2 with 
the three solutions 

f2 = iO,±^V3XS^). (7) 

'We remark that, in addition to the compacton solutions discussed below, Eq. (|5]l has the iso- 
lated (i.e., without free integration constants) analytical solution (p ~ ±y/X/2Qs^. This solution 
grows without bound for increasing s and, obviously, cannot give rise to a finite energy configura- 
tion. 
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If we further assume, without loss of generality, that takes the positive value 
/2 = +(3/2)\/3AS"^ for s < 5 and the vacuum value /2 = for s > S then the 
higher coefficients fj are uniquely determined by linear equations and we find the 
expansion 

= -y/3XS\s-Sf + —y/3XS^{s-Sf + 
2 5 



2V3XS\s - S)^ + 0{{s - Sf) for s<S 
(f){s) = for s>S. (8) 

We remark that also in the r variable the function 0(r) approaches its vacuum 
value quadratically at r = i? = S^, because the conditions (p{S) = (p'{S) = are 
invariant under a variable change. Indeed, we find 

r-R)^ + 0{ir-R)^), for r<R 



270i?2 V 3 
0(r) = for r > R. 

2.2 Expansion at the center 

Plugging a Taylor series expansion around s = into eq. (|5]) we find that the 
first two coefficients </)(s = 0) = 0o and 0'(s = 0) = a remain undetermined, 
whereas the higher ones are determined uniquely by linear equations. Due to the 
strong suppression factor at the r.h.s. of eq. Q we find that the coefficients of 
to are, in fact, zero. Concretely we find 

= h + as + ^A.i° + ^Xs'' + (9) 

The three next higher order terms 0{s^^) — 0{s'^^) include only terms proportional 
toA^. 



2.3 Compact ball solutions 

We shall find solutions explicitly by numerical integration. But before performing 
these calculations we want to discuss some generic features of eq. ([5]) which allow 
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to establish both the existence and several properties of these compact ball solu- 
tions. A first indication for the existence of solutions can be found by counting the 
number of boundary conditions and free parameters in the numerical integration. 
There are two possibilities for the numerical integration, namely a shooting from 
the center, or a shooting from the boundary. If we shoot from the center, there 
are altogether three free parameters, namely 0(0) = 0o, = ^'^d the com- 
pacton radius 5*. On the other hand, there are two boundary conditions that have to 
be obeyed for a compacton, namely 0(5*) = 0'(5') = 0. Therefore, we expect the 
existence of a one-parameter family of compact balls which may be parametrized, 
e.g., by the compacton radius S. If we shoot from the compacton boundary, the 
only free parameter is the compacton radius S. On the other hand, there are no 
boundary conditions at s = 0, because both 0(0) and 0'(O) are undetermined by 
the expansion at s = 0. So again we expect a one-parameter family of solutions. 

Further conclusions may be drawn by a closer inspection of eq. ([5]). Firstly, 
we observe that eq. ([5]) is invariant under the reflection —0. Therefore 
we may choose 0o = 0(0) > without loss of generality. Secondly, it follows 
immediately from eq. ([5]) that away from the center (i.e. for s > 0) it holds that 
if > then 0" > 0. As a consequence, if 0o > and a = 0'(O) > then 
the scalar field will grow indefinitely and may never reach its vacuum value 
= 0. We conclude that a compact ball requires a < 0. Now for generic a < 
the following may happen. If for a fixed value of 0o the slope is too weak (i.e. 
if |a| is too small) then 0(s) will reach a point sq where 0'(so) = before it 
reaches = 0. At this point sq the second derivative 0" (and, consequently, all 
higher derivatives) becomes singular and the integration breaks down. On the 
other hand, if |a| is too big, then will cross the line = instead of touching 
it, and then develop towards more and more negative values of indefinitely. It 
follows that there exists an intermediate value of a such that touches the line 
= instead of crossing it. This is the compact ball solution, and the point s = S 
where the touching occurs is the compacton radius. 

The explicit numerical integration confirms the behaviour described above to 
a high precision. In the numerical integration we mainly used shooting from 
the compacton boundary because of the lower number of free parameters, which 
makes the scan for compacton solutions less time consuming. A typical solution 
of the numerical integration is shown in Figs. 1, 2. Here, in Fig. 2 we plot the 
solution of Fig. 1, but using r as the independent variable. The non-analytical 
behaviour in r (the spike at r = 0) can be clearly seen. 

Finally, let us remark that, as always in the case of compactons, we may ob- 
tain multi compacton configurations by distributing several non-overlapping com- 
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pactons with different centers in flat space. 



2.4 Energy functional 

We found that there exist compact ball solutions for arbitrary values of the field 
at the center 0o (oA equivalently, for arbitrary values S of the compacton radius), 
and it is easy to see that these different values of 00 correspond to different values 
of the total energy. Therefore, compact balls do not correspond to genuine criti- 
cal points of the energy functional, and one may wonder why solutions do exist 
at all. There are two possible ways to understand this puzzle. Starting from the 
reduced energy functional for radially symmetric fields, the answer is that varia- 
tions of the energy functional may receive contributions from the boundary, and it 
is precisely these variations of the boundary which give nonzero contributions to 
the energy functional for a compact ball solution. The second, equivalent answer 
is that the compact balls solve an inhomogeneous equation with a delta function 
source term. This we want to demonstrate explicitly in the sequel. In a first step, 
we write down the energy functional and then use the principle of symmetric crit- 
icality to rewrite this energy functional as a functional for radially symmetric field 
configurations (the principle of symmetric criticality just states that this reduction 
of the energy functional to radially symmetric configurations provides the correct 
radially symmetric field equations if the field equations are compatible at all with 
this symmetry reduction - i.e., with the ansatz (f) — (t){r)). We get 



E[(t)\ = J r'^dr sin eded^(^^{(y(pYy + X(j)^ 

= An l\^dr(^^(l)^^ + X(l)^y (10) 



For the variation of the functional we find 

5E = An I r^dr ^ 







r 



= An p dr (-dr{r'^(t)l) + 2r^A0) (50 

+An{r^(ty'M)\o (11) 

where we performed a partial integration and used the fact that is a compact ball 
solution which takes its vacuum value for r > R. Using the small r behaviour 
0(r) ~ 00 + cer 3 we finally get for the boundary term 

47r(rV'50)|^=-47r— 50(0). (12) 
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Therefore, the variation of the energy functional for compact balls is zero only 
under variations which do not change the value of at the center, 5(f)(0) = 0. 

A second, equivalent interpretation is like follows. We may require invari- 
ance of the energy functional under general variations, but then, using 5(f)(0) = 
J dr5(r)5(f)(r), we find the inhomogeneous equation with a delta function source 
term, 

- dr(r^(Pl) + 2r2A0 - —5(r) = 0. (13) 

Here the strength of the source term is related to the value of at the origin, 0o, so 
different 0o (and, therefore, different energies) just correspond to different equa- 
tions, i.e., different strengths of the source term. The coefficient of the source 
term depends directly on a = (f)'(0), but for finite energy solutions (that is, com- 
pact balls), there exists a definite relation between a and 0o, so the source term 
determines 00 uniquely. 




10 r 




Figure 1: Compacton without gravity - the profile of the scalar field (f>(s) and its 
derivative 0'(s) for A = 1. Compacton radius S = 1.5. Values at the center 
0(0) = 8.271, 0'(O) = -6.355. 
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Figure 2: Compacton without gravity - the profile of the scalar field of Figure 1, 
but expressed in the variable r = s^, 0(r). Compacton radius R = = 3.375. 



3 The model with gravity 

We now consider the scalar field model of the previous section coupled minimally 
to gravity. The action reads 

S = J d^x^l (^7^ -X\X\- 1^(0)) , (14) 

where, as before, 

X = \g^''d,4>d,(P, V{4>) = x<p\ (15) 

Further, IZ is the curvature scalar, and g is the determinant of the metric tensor. 
We now assume a spherically symmetric space time, then the line element ds^ 
may be chosen in the form 

ds"^ = -A{r)dt^ + B{r)dr'^ + r\de'^ + sin^ ed(j)'^) (16) 

where for consistency we have to assume that (p = (pir) depends on the radial 
coordinate r only. The Euler-Lagrange equation for 0(r) takes the form 

+ + ^ + f )) \B-'\B-^^'' - 2A0 = 0. (17) 
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The Einstein equations may be obtained from the Einstein tensor, which has the 
independent components 

Goo = A[irBY'B' + ^il-B-')y (18) 



Gu = B(^irAB)-'A'-^{l-B-')y (19) 
G22 = (^^iAB)-'/\iAB)-'/^Ay + ^irAB)-'A' -^{rBY'B'yiO) 

(G33 is not independent but obeys G33 = sin^ OG22), and the energy-momentum 
tensor 

dC 

= 2\X\d^<pd,(P-g^,i\X\X + X(j)^) (21) 
with the independent components 

Too = aQ|5-1|5-V" + A02), (22) 
Tn = \B-^\4>"'-b(\\B-^\B-^^'^ + X^A, (23) 



A 

T22 = -r2('i|5-i|fi-V" + A02). (24) 



The Einstein equations read 



B' = -B{1-B) + k\B^(-\B-^\B-^(I)'^ + X(I)^), (25) 



r V4 

^ = -(B-l) + -K^r\B-^\(j)'^-K^XrB(j)\ (26) 
A r 4 

^{AB)-^A" - ^{AB)-\AByA' + ^{rAB)-^A' - ^{rB^y^B' + 

+«:2Q|fi-i|B-V' + A02) =0. (27) 

Here several comments are appropriate. Firstly, there seem to be four equations 
(the field equation (flTI) and the three Einstein equations) for three real functions 
(p, A and B. As always in the case of one real scalar field, these equations are, 
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however, not independent. The field equation may, in fact, be derived from the 
three Einstein equations. Further, in the case at hand, the field equation (fTTl ) may 
be derived from the three Einstein equations in a purely algebraic fashion (i.e., 
without performing additional derivatives). This follows easily from the fact that 
both the field equation (flTl) and the third Einstein equation (ITTl) are of second 
order. Therefore, the field equation and the third Einstein equation are completely 
equivalent, at least in non-vacuum regions where 0' 7^ (in regions where 0' = 
the field equation is more restrictive, because its derivation from the Einstein 
equations requires a division by 0')- We remark that the above feature is not 
always true (in some cases the Einstein equations are more restrictive and have to 
be used, because the derivation of the field equations involves further derivatives). 

Secondly, the function A appears in all equations only in the combination 
{A' /A) and derivatives thereof (this is also true for Eq. (l27l) . where it is not com- 
pletely obvious). Therefore, one may eliminate the function A with the help of Eq. 
(I26l) from all the remaining equations. One may choose a set of two independent 
equations in (p and B from these remaining equations, solve them for (p and B, 
and determine the corresponding A from Eq. (|26l) in a second step. It is obvious 
from Eq. (|26l ) that A is determined up to a multiplicative constant. The choice of 
this constant just corresponds to a constant rescaling of the time coordinate. This 
constant will be fixed by the condition that asymptotically (that is, for r bigger 
than the compacton radius R ), the metric is equal to the Schwarzschild metric, 
i.e., 

A{r) = B-\r) = 1 - — for r > R > Rs (28) 

r 

(we will find that the compacton radius is always bigger than the Schwarzschild 
radius, R > Rg, and that a horizon never forms). Concretely, for the two indepen- 
dent equations for B and (p we choose Eq. (|25l) for B and the field equation Eq. 
(flTI) for (f), where we eliminate both (A' /A) and B' from the latter equation with 
the help of Eqs. (|26l) and (l25l) . It will be useful for our discussion to display the 
resulting system of two equations again. We get 

f = + + (29) 

for B and 

= 2 (^K^Xr(f)^B - 0'^ + 2Asign(5)5V (30) 

for (p (here sign(i?) = (\B\/B) is the sign function). We remark for later use that 
if we forget about the sign function, then the equations for negative B may be 
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recovered from the equations for positive B by the combined coupling constant 
transformations 

A^-A, K^-*-K^. (31) 
3.1 Qualitative behaviour of B 

Before starting the numerical investigation and the expansions at the compacton 
boundary and at the center, we want to draw some conclusions on the behaviour 
of the function B. Concretely, we want to show that for a nonsingular scalar field 
0, B(r) cannot approach zero from the inside, that is, from smaller values of r. 
This implies that when B{r) takes the value zero at some radius r = tq, then 
B{r) is not defined for r < tq. Differently stated, if we start the integration at 
some value r > tq where B{r) > (this we assume because we want to connect 
to the Schwarzschild solution), and we then integrate downwards (i.e. towards 
smaller values of r), then the integration breaks down at r = tq. We shall find in 
the sequel that at this value a singularity forms, that is, some curvature invariants 
become infinite at r = tq. Here we have to distinguish two cases, namely rg = 
or To > 0. In the case tq = the singularity is just a point at the origin of our 
coordinate system, whereas for rg > the locus of the singularity is a two-sphere 
5"^, so space has a singular boundary with S"^ topology. We shall call solutions 
of the first type "compact balls" in the sequel, whereas the second type (with the 
singular inner boundary of shape) are called "singular shells". We discuss this 
behaviour here because it is precisely what is found in the numerical integration 
(that is, if we start at some compacton boundary r = R and then integrate towards 
smaller r, B will always hit the line B = 0, either at r = or at some nonzero 
r = To). 

It remains to demonstrate that B cannot be continued to the region r < tq. If 
-B > for r > To and approaches zero at r = tq, then necessarily {B' / B) > 
for r > ro but sufficiently close to tq. A continuation to r < tq may either cross 
zero, in which case B < and 5' > for r < ro but sufficiently close to tq. 
Or the continuation may return to the region B > for r < ro, in which case 
B > and B' < for r < tq but sufficiently close to tq. In both cases the ratio 
between B' and B is negative, {B' / B) < for r < rg but sufficiently close to ro. 
The inequality {B' / B) < is, however, incompatible with Eq. (|29l ) for r < ro 
and sufficiently close to ro. In fact, the r.h.s. of Eq. (|29l ) is manifestly positive 
for r close to ro. The first term (l/r)(l — 5) ~ (1/?") is positive and not small. 
The second term, K^r{(j)'^ / A\B\), is positive, and its magnitude depends on (p. The 
third term, K^XrBcj?, may change sign if we assume that B does. It is, however, 
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small near tq, because it is proportional to B itself. The conclusion is that the 
ratio (B' / B) is always positive near r = tq and, therefore, B cannot be continued 
to values r < tq. 

In a next step, we want to prove that a Schwarzschild type horizon of B is im- 
possible for compacton solutions, too. We remind the reader that the Schwarzschild 
solution for B is 

Bs=(l--y' = ^^ = - + 1 , e^r-R,. (32) 
V r / r — Rs e 

For a Schwarzschild type horizon coordinate singularity we, therefore, expand 

p oo 

^ n=0 

oo 

</'(e) = E/n^" (33) 

n=0 

We now insert this expansion into Eq. (l30l) where we assume at the moment that 
r > Rs, that is, e > and, therefore, B > 0. Eq. (|30l) then becomes 

= 2 f fi:'A(i?, + e)0'5 ^) 0'^ + 2A52(/, (34) 

V Rs + eJ 

and we immediately find that /o = because otherwise the coefficient of e^^ in 
Eq. (|34l) cannot be set to zero. For /i we find a cubic equation with the three 
solutions 

/i = (0,±i?,v^) (35) 

For /i = it follows easily that all the higher fi are zero, as well, so that = 
takes its vacuum value everywhere, and we are back to the Schwarzschild solution 
B = Bg for the metric function B. For non-zero /i we may choose the positive 
root /i = RsVX without loss of generality. This choice corresponds to a formal 
solution, but the resulting scalar field grows without bound for increasing e, as 
follows easily from Eq. (I34l) . Indeed, can avoid unbound growth only if 0' = 
somewhere in the region e > 0, but 0' = implies 2A-B^0 = which is impos- 
sible (0 is nonzero at its maximum by assumption, whereas B cannot approach 
zero in the direction of growing r, as we know from the last paragraph). Specif- 
ically, the formal solution for nonzero /i can, therefore, never join a compacton 
with its compacton boundary somewhere in the region r > Rg. Reversing the 
argument, we conclude that a horizon of the Schwarzschild type can never form 
for a compact ball solution. 
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We still want to know what happens near a Schwarzschild type horizon for 
r < Rs, i.e., e < 0. There B < 0, and we remind the fact that solutions for _B < 
may be inferred from solutions for i? > by the coupling constant transformation 
(|3TI) . specifically A — *• —A. We therefore find for the linear coefficient /i = 
(0, ±iRs^JX), and obviously only the trivial vacuum solution /i = =^ = 
is real and, therefore, physically acceptable. We conclude that the only possible 
field configuration inside a horizon is the vacuum configuration </> = 0. We might 
still ask what happens if we try to put a compacton strictly inside the horizon, that 
is, at a compacton radius R < Rg. We will find that the result is the same, i.e., 
the nontrivial expansion coefficients for become imaginary, and only the trivial 
vacuum configuration = is possible inside the horizon. 

We conclude that for compact boson star solutions the function B can never 
be negative, B > 0. 

3.2 Behaviour at the boundary 

Now we assume the existence of a compact boson star boundary, that is, a radius 
r = R where 0(-R) = 0'(-R) = 0, whereas the second derivative is zero from 
above but nonzero from below. Assuming that B is nonnegative (because we 
want to smoothly join it to the Schwarzschild solution for r > R), and plugging 
the power series expansions around r = R 



(36) 



k=2 



R)\ 



(37) 



into (Ull) and dSO]) gives 




h 



'1 



R 



{bo - 1) 
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b2 



As expected, we find three roots for /2, so the vacuum solution /2 = for r > R 
may be smoothly joined to one of the two nontrivial roots of /2 at r = R. In the 



sequel we choose the positive root /2 = (l/2)y A/36o without loss of generality. 
We remark that for negative B the nonzero roots of the coefficient /2 become 
imaginary, as already announced in the previous section. 

Inserting, further, the power series expansions above and the one for A, 



into Eq. (|26l) for A we get 



A{r) = ^ak{r-Rf (38) 

k=0 



For the expansion coefficients of B and A we, therefore, find that also for nonzero 
fk, the leading behaviour of A{r) and B{r) close to r = is like for the vacuum 
(i.e., Schwarzschild) solution. The first contribution of a nonzero /2 appears in ag 
and in 65. An easy way to see this is to study the power series expansion of the 
product AB for r < R, 

Air)B{r) = ao&o - ^n'X^aob^iR - r)' + OUR - r)^), (39) 

where the positive nontrivial root of /2 was inserted. For r > R, A and B should 
form the Schwarzschild metric, which fixes the constant oq to the value = 69 ^. 
On the other hand, bo is a free parameter which is related to the Schwarzschild ra- 
dius Rg or to the asymptotic Schwarzschild mass = {Rs/'^) of the asymptotic 
Schwarzschild metric, as well as to the compacton radius, via 
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We find that the expansion at the compacton boundary leaves us with two free 
parameters, namely and the compacton radius R. Equivalently, we may choose 
the Schwarzschild mass rus = Rs/2 and the compacton radius as free parameters. 

Before studying the expansion from the inside, let us discuss briefly what to 
expect for an integration which starts at the compacton boundary and proceeds 
towards smaller r. It holds that at r = R, B{R) > 1 and (j){R) = (p'i^R) = 0, 
therefore B'{R) is negative and B will at first increase towards smaller values of 
r, see Eq. (|29l ). For even smaller values of r, the additional, positive terms at 
the r.h.s. of Eq. (|29l ) start to contribute, so B' may become positive or negative, 
depending on the relative strength of the different terms at the r.h.s. of Eq. (|29l ). 
It turns out numerically that for a certain radius B' always becomes positive such 
that B starts to shrink towards smaller values of r. Once B has shrunk sufficiently 
such that B < \, then the r.h.s. of Eq. (|29l) is necessarily positive, such that B has 
no other choice than shrinking further. Numerically, this is exactly what happens, 
for all possible values of the two free parameters. Therefore, for sufficiently small 
r, there are the following three possibilities for the behaviour of B. It may shrink 
to a nonzero value at r = 0, i.e., B{r = 0) > 0. We will see that there exists only 
one isolated solution for 5(0) > (i.e., without free integration constants), and 
this solution cannot be connected to a compacton boundary. Or B may go to zero 
atr = 0,-B(r = 0) = 0. We will see that two different kinds of solutions of this 
type may be connected to a compacton boundary. They will form the solutions of 
the compact ball type. The third possibility is that B becomes zero already for a 
nonzero r = tq, i.e., B{r = ro > 0) = 0. These are the singular shells. 

3.3 Expansion at the center 

In this section we assume that a solution exists locally near r = 0. We insert the 
power series expansions for and B 

0(^) = T.hr\ (41) 

fe=0 

B{r) = Y.bkr\ (42) 

k=0 

into (|29l ) and (|30l ). In a first step we want to assume that 6o > 0. Cancellation of 
the coefficient of r^^ in Eq. (|29l ) then requires Iiq = 1, whereas cancellation of 
the coefficients of r^^ and r*^ in Eq. (|30l) require /o = and /i = 0. All higher 
expansion coefficients are uniquely determined (up to an overall sign of 0), that 
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is, the solution is an isolated one without free integration constants. Explicitly, we 
find (we choose the plus sign for 0) 



<p(r) = ^r^ + 0(r^ 
a/20 



B{s) = ^+^-^r' + 0{r'') (43) 

This solution behaves like at a compacton boundary already at the center r = 
(i.e. 0(0) = 0, 0'(O) = 0), and the question is whether it can be connected to 
a compacton boundary at some nonzero r = _R. It follows easily from Eq. (|30l) 
that this is impossible. Indeed, if takes its vacuum value = at two different 
radii, then it must pass through a local maximum at some r = tq between the two 
radii (where 0(ro) > and 0'(ro) = 0). But the existence of this maximum is 
incompatible with Eq. (|30l) . 

Therefore, we now assume 6o = and find 



B{s) = hr+{]K^f^-bi)r^+(bl-^K%f^)r^ + 0{r^). (44) 



3 V27 36 

Here, /q, /i and bi are free parameters (integration constants). We remark that for 
compact ball solutions (that is, for the correct matching to a compacton boundary) 
both /o and /i have to be nonzero. On the other hand, there will exist solutions 
with bi = 0. In this case the leading terms of the expansion read 



jo^ji ggji -r 2160 ^ 



- r -t- L^i 

B(s) = ^!/y_!^!S.^ + f!fill±M).e + 0(r') (45) 

^ ' 4 •^i 144 432 K J \ J 



The power series expansion for A has to be treated independently for bi ^ and 
bi = 0, because the leading behaviour for small r is completely different, being 
A(r) ~ r^^ for bi ^ 0, and A ~ for bi = 0, as may be inferred easily from Eq. 
(I26l) . For bi ^ the leading terms in the expansion read 

Air) = a., (i + U + ^) + ^ (46? + 3.^7) r + 0(r^)) (46) 
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whereas for 61 = it reads 

A{r) = «2 (r^ - + + 0(r«)) (47) 

The leading coefficients a_i or 02, respectively, are free parameters in the power 
series expansion. They appear as linear factors also in all higher coefficients, in 
accordance with our observation that A is determined up to a multiplicative con- 
stant. They cannot be determined from the local analysis at r = 0, but must instead 
be determined from the condition that A approaches a Schwarzschild metric at the 
compacton boundary. 

3.4 Selfgravitating compact ball solutions 

Before explicitly performing the numerical integration, we again want to deter- 
mine the number of free parameters and boundary conditions, both for a shooting 
from the boundary and for a shooting from the center. Here we only consider the 
parameters and boundary conditions for (p and B, because these two can be deter- 
mined from the system of two equations (|29l ) and (|30l) . A may then be determined 
from Eq. (|26l) in a second step, and we know already that the multiplicative con- 
stant which is the free parameter of A must be determined by a matching to the 
Schwarzschild metric at the compacton boundary. 

At the compacton boundary there are two free parameters, namely the com- 
pacton radius R and the Schwarzschild (or asymptotic, ADM) mass m^. Con- 
cerning the conditions that have to be imposed at the center r = 0, we have to 
distinguish the case bi ^ from the case 61 = 0. In the case 61 7^ 0, there are 
no conditions imposed at r = 0, because the first two coefficients /o and /i of (f) 
are unrestricted, and the fact that 60 = does not count as a boundary condition, 
because it is true for all generic solutions that exist locally near r = (we remind 
the reader that the solution with 60 = 1 is an isolated solution with no free param- 
eters). We, therefore, expect a two-parameter family of solutions. Specifically, for 
a fixed Schwarzschild mass we expect to find a one-parameter family of compact 
ball solutions with different radii. We shall call this type of solutions with 61 > 
"large compactons" in the sequel. 

In the case bi = 0, this condition provides exactly one boundary condition 
at r = 0. In this case we, therefore, expect a one-parameter family of solutions. 
Specifically, for a fixed Schwarzschild mass rus we expect only one compact ball 
with a fixed radius. We shall refer to this type of solutions as "small compactons". 
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An analysis of the shooting from the center leads to the same results. For the 
case bi ^ (large compactons), there are four free parameters, namely bi itself, 
/o, /i, and the compacton radius R. Further, there are two boundary conditions at 
the compacton boundary r = R, namely the conditions = and (j)'{R) = 0. 
Therefore, we expect a two-parameter family of compact ball solutions. In the 
case bi — (small compactons), we are left with three free parameters /o, /i 
and R and the same boundary conditions atr — R, therefore we expect a one- 
parameter family of solutions. 

The explicit numerical integration completely confirms the above results. Con- 
cretely, we prefer to shoot from the boundary because of the smaller number of 
free parameters. There, the free parameters are the Schwarzschild mass = 
Rs/2 (where Rg is the Schwarzschild radius), and the compacton radius R. In the 
numerical integration we choose a fixed and then perform the integration for 
different values of R. We already know from the discussion of the previous sec- 
tions that a solution cannot exist for R < Rs. For i? > we find the following 
behaviour. If R is too small, then in the integration from r — R towards the center 
B reaches zero already at a nonzero radius r = rg, i.e., i?(ro) = 0. These singular 
shell solutions shall be discussed in the next subsection. Increasing R further, it 
reaches a minimum value R = Rgc (where sc stands for "small compacton") such 
that B approaches zero at r = 0. For this minimum value of R, B approaches zero 
quadratically, B ~ r^, i.e., it corresponds to the case bi — 0, the small compacton. 
Two examples (for nis — 0.1 and = 1.5, respectively), are shown in Figures 
3-8. 

Further, we show, in Figures 9-11 the dependence of R and the values of 
and (f)' at the origin r = as a function of the Schwarzschild mass rus (which 
we use as the only free parameter in the case of small compactons). We remark 
that for sufficiently large Schwarzschild mass m^, the relation between the com- 
pacton radius Rgc and the Schwarzschild mass plotted in Fig. 9 may be well 
approximated by the linear interpolation formula 

Rsc ~ 2.9ms + 3.4 > 2ms, (48) 

and Rsc is always above the Schwarzschild radius Rg = 2ms, i-e., Rsc > Rs 
(this inequality remains true also for small values of ms where the above linear 
interpolation cannot be used). 

For values of R which are larger than Rg^ for a given mg (the radius of the 
small compacton), B approaches zero linearly. This corresponds to the case bi > 
0, i.e., to the large compacton. There exist solutions of this type with arbitrarily 
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Figure 3: Compacton with gravity included, small compacton case (bi = 0). 
Profile of the scalar field and its derivative for A = 1 and k = 0.1. Shooting from 
the boundary: initial value = 0.1, resulting compacton radius Rgc = 3.415. 
Values at the center 0(0) = 7.633 and 0'(O) = -14.341. 

B 

1.5 r 




L 0.5 1.0 1.5 2.0 2.5 3.0 

Figure 4: Small compacton case (bi = 0). Function B{r) for A = 1 and k = 0.1. 
Shooting from the boundary, like in Figure 3: = 0.1, Rsc = 3.415. 



large compacton radius R for one and the same value of the Schwarzschild mass. 
Some examples are shown in Figures 12-14 for the fixed value = 1.5 and three 
different values of the compacton radius. It can be seen from Figures 13, 14, that 
the functions A and B approach the same Schwarzschild solution (with the same 
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Figure 5: Small compacton case {hi = 0). Function A{r) for A = 1 and k = 0.1. 
Shooting from the boundary, like in Figure 3: nis = 0.1, Rgc = 3.415. 








Figure 6: Compacton with gravity included, small compacton case (bi = 0). 
Profile of the scalar field and its derivative for A = 1 and k = 0.1. Shooting from 
the boundary: initial value = 1.5, resulting compacton radius Rgc = 7.800. 
Values at the center 0(0) = 12.021 and 0'(O) = -1.481. 



Schwarzschild mass = 1.5) for sufficiently large radius. There exist, in fact, 
some subtleties related to the large compacton solutions, which we want to discuss 
now. First of all, smaller values of bi do not correspond to smaller values of the 
compacton radius (specifically, the radius of the small compacton is not reached 
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Figure 7: Small compacton case (bi = 0). Function B{r) for A = 1 and k = 0.1. 
Shooting from the boundary, like in Figure 6: rus = 1.5, Rsc = 7.800. 
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Figure 8: Small compacton case (bi = 0). Function A (r) for A = 1 and k = 0.1. 
Shooting from the boundary, like in Figure 6: = 1.5, Rsc = 7.800. 



by the limit 6i ^ of the large compacton radius). Quite on the contrary, small 
bi corresponds to large compacton radius, as can be seen easily in Figure 13. This 
demonstrates again that the small compacton case bi = cannot be seen as the 
limiting case lim 6i ^ of the large compacton case. 

A second subtlety can be seen once we reverse the role of nis and R in the case 
of the large compactons. That is to say, we choose now a fixed R and then vary 



22 



0.5 



1.0 



1.5 



2.0 



Figure 9: Small compacton case (bi = 0). Radius Rgc of the small compacton as 
a function of the Schwarzschild mass m,. 
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Figure 10: Small compacton case (bi = 0). Value of the scalar field at at the 
center, as a function of the Schwarzschild mass, ff^inia). 



nis- For too large m^, no solution exists. For a certain limiting value of m^, we get 
the small compacton, and for even smaller m^, we find the large compactons. The 
question arises how small mg can get for a given R such that a large compacton 
exists. The answer is that it can get arbitrarily small. Specifically, there exist 
large compactons for negative values of the Schwarzschild mass m^. We show an 
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Figure 11: Small compacton case (bi = 0). Value of the first derivative of the 
scalar field at at the center, as a function of the Schwarzschild mass, f^{ms). 




Figure 12: Compacton with gravity included, large compacton case (bi > 0). 
Profile of the scalar field and its derivative for A = 1 and n = 0.1. Shooting 
from the boundary: initial value = 1.5, and the three initial values for the 
compacton radius R = 9.3, R = 10.3 and R = 12.8. 



example of this fact in Figures 15-17. 

The existence of large compacton solutions for a negative Schwarzschild mass 
may be further understood by the introduction of a "variable mass" function m(r) 



24 



B 

2.0 r 




Figure 13: Large compacton case (61 > 0). Function B{r) for A = 1 and k = 0.1. 
Shooting from the boundary, like in Figure 12: = 1.5, and R = 9.3, R = 10.3 
and R = 12.8. 




Figure 14: Large compacton case (61 > 0). Function A{r) for A = 1 and k = 0.1. 
Shooting from the boundary, like in Figure 12: = 1.5, and R = 9.3, R = 10.3 
and R = 12.8. 



as a new variable instead of B, 



B(r) = 1 



2m(r) 



-1 



r — 2m(r) 



(49) 
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Figure 15: Large compacton case (61 > 0). Profile of the scalar field and its 
derivative for A = 1 and k = 0.1. Shooting from the boundary for the initial value 
R = 6 and the three initial values rus = 0.5 (upper curve), = (middle curve) 
and nis = —0.5 (lower curve). 
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Figure 16: Large compacton case (61 > 0). Function B{r) for A = 1 and k = 0.1. 
Shooting from the boundary for the initial value R = 6 and the three initial values 
nis = 0.5, = and m^, = —0.5. The graph for positive has a maximum 
and decreases for larger r, the graph for = approaches the constant value 1 
at r = i?, whereas the graph for negative never reaches B = 1 and has positive 
slope for all r. 
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Figure 17: Large compacton case (hi > 0). Function A{r) for A = 1 and k = 0.1. 
Shooting from the boundary for the initial value R = 6 and the three initial values 
nig = 0.5 (lower curve), iris = (middle curve), and nis = —0.5 (upper curve). 



The condition that B{r = 0) = just requires that m(r) should not go to zero 
too fast in the limit r ^ (m(r = 0) 7^ is a sufficient condition), whereas the 
condition B'{r = 0) =bi > leads to 

mo = m{r = 0) = — ^ < 0. (50) 

This relation will explain the unexpected relation between bi and the radius R 
of the large compacton. In addition, this relation demonstrates, once again, that 
the small compacton case cannot be obtained as the limit of the large compacton 
case (the function m(r) is, in fact, not useful for the small compacton, because it 
does not have a power series expansion at r = for the small compacton case). 
Inserting m(r) into Eq. (|29l ) for B we get 

m' = y ((r - 2m) + Ar^^^) (51) 

with the formal solution 

m(r) = / df — {{f — 2m)^0'^ + Af^0^ ) + uiq = m{r) + mo (52) 
JO 2 ^ ^ 

and, therefore, for the Schwarzschild mass 

nis = m{r = R) = m{R) + mo. (53) 
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We observe that the integrand in the definition of m(r) is non-singular and pos- 
itive definite, therefore m{R) is a monotonously increasing function of R and 
obeys rh(R = 0) = 0. Now we treat R and mo (that is, R and bi) as our two 
free parameters for the large compacton solutions. The only restrictions on these 
parameters are R > and mo < (that is, bi > 0). Obviously, for any choice of 
R > there always exists a choice for mo < such that the Schwarzschild mass 
is negative, m^ < 0. One may, of course, exclude such choices as unphysical. We 
just want to remark that from the point of view of the large compacton solutions 
there is nothing special about negative m^, i.e., solutions with negative nig may 
be reached from positive m^ by a completely smooth variation of the parameters 
which characterize these solutions. 



3.5 Singular shells 

Now we want to discuss the case of singular shells, that is, of solutions which 
connect a singular inner boundary at r = tq with a compacton boundary atr = R. 
Here, the singular boundary is defined by B(ro) = 0, and it holds that ro < R. In 
a first step, we insert the power series expansions about r = ro 

0W = T.M^-^o)\ (54) 

k=0 

B{r) = Y.Ur-To)\ (55) 

k=l 

into (|29] ) and (l30l) . where we assume B{r = ro) = 0. We get 

= /o + /i(r-ro) + ^(«:2Aro7^-l)(r-rof + O((r-ror) 
B{s) = l«:Vo/7(r-ro) + ^«:7i(r-ro)2 + 

^«:Vo/7 (n'Xrlf^) (r - ro)^ + 0((r - ro)^). (56) 

Here, /o and /i are free parameters (integration constants), whereas bi is already 
determined by the equations. Inserting the corresponding power series expansion 
for A into Eq. (|26l ). we find that the first nonzero coefficient is the cubic one in 
this case. Explicitly, we get 

A{r) = a^fi^r - rof - —{r - roY+ 
V ro 

+ ^ {^'>^rlf^ - l)) - ^o)' + Oiir - ro)«)) . (57) 
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Next, let us count free parameters and conditions in this case. For a shooting from 
the outer compacton boundary, we have three free parameters, namely m^, R and 
tq. On the other hand, we have one condition, namely B(ro) = 0. So we expect a 
two-parameter family of solutions. If we shoot from the inner boundary, we have 
the four free parameters /o, /i, R and tq, and the two conditions (f){R) = (p'{R) = 
0. Again, we expect a two-parameter family of solutions. Specifically, for a fixed 
asymptotic mass we expect a one-parameter family of solutions which may be 
parametrized, e.g., by the compacton radius R. This is in complete accordance 
with the discussion of the previous section where we found already that for fixed 
Rs = 2ms there should exist singular shell solutions for Rg < R < Rgc (here Rsc 
is the radius of the small compacton). It turns out that the inner singular radius of 
the singular shell is always smaller than the Schwarzschild radius, so we have, in 
fact, tq < Rs < R for all singular shells. 



10 - 




Figure 18: Singular shell case. Functions 0(r) and 0'(r) for A = 1 and k = 0.1. 
Shooting from the compacton boundary for the initial values = 0.5 and R = A. 
The resulting singular boundary is at tq = 0.4736 (position of the minimum of 

0'). 

We show two cases in Figures 18 - 25, where we choose Rg = 2ms = 1 in 
both cases. The corresponding small compacton radius is Rsc = 4.885, and we 
choose _R = 4 in Figs. 17-21, and i? = 2 in Figs. 22 - 25. For a value of R near 
Rsc (the first case R = 4), the metric functions deviate form the Schwarzschild 
behaviour rather soon, and the transition from r = i? to r = ro is rather mild. 
For a R near Rg, on the other hand, the metric functions B and A behave almost 
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Figure 19: Singular shell case. Function B{r-) for A = 1 and n = 0.1. Shooting 
from the compacton boundary for the initial values nis = 0.5 and R = A. The 
resulting singular boundary is at tq = 0.4736 (position of the zero of B). 
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Figure 20: Singular shell case. Function A{r) for A = 1 and k = 0.1. Shooting 
from the compacton boundary for the initial values = 0.5 and i? = 4. The 
resulting singular boundary is at tq = 0.4736 (position of the zero of A). 



exactly like the Schwarzschild solution for most values of r. This is clearly seen 
in the plot of the product AB in Figure 25. Specifically, B almost reaches the 
Schwarzschild horizon divergence (in the concrete example for i? = 2, it has a 
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Figure 21: Singular shell case. Product A{r)B{r) for A = 1 and k = 0.1. Shoot- 
ing from the compacton boundary for the initial values = 0.5 and i? = 4. The 
resulting singular boundary is at tq = 0.4736 (position of the zero of AB). 



maximum value of -Bmax > 500), before steeply descending to zero, see Fig. 23. 

We shall find in the next section that the condition that B{ro) = at the inner 
shell boundary r = tq indeed implies that some curvature invariants become sin- 
gular there, justifying the name "singular shell". In fact, space-time itself ends at 
the singular boundary, i.e., cannot be continued beyond the singularity (remember 
that we found in Section 3.1 that there is no continuation of B(r) to values r < tq 
if -B(ro) = 0). We will also find, in the next section, that the singularity is, at 
the same time, the locus of a Killing horizon, and no signals can propagate from 
the singular surface to an outside observer, i.e., the singularity is invisible. In this 
sense, these singular shells have some similarity with ordinary black holes. 

3.6 Singularities of curvature invariants and radial geodesies 

We found two types of solutions in the last sections, namely compact balls and 
singular shells. Within the compact balls, solutions are further discriminated by 
the approach of B and A to the center r = 0. The first case of large compactons 
is the case bi > 0, where B approaches the center linearly, B ~ bir, and A ~ 
a_ir~^. The second case of small compactons is defined by 6i = and, as a 
consequence, B ~ 62^"^ and A ~ a2r^ (here 62 is, in fact, not a free parameter, but 
this issue is irrelevant for the present discussion). 
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Figure 22: Singular shell case. Functions 0(r) and 0'(r) for A = 1 and k = 0.1. 
Shooting from the compacton boundary for the initial values = 0.5 and R = 2. 
The resulting singular boundary is at tq = 0.9928 (position of the minimum of 0', 
which is outside the range of the plot). 
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Figure 23: Singular shell case. Function B(r) for A = 1 and k = 0.1. Shooting 
from the compacton boundary for the initial values = 0.5 and R = 2. The 
resulting singular boundary is at tq = 0.9928 (position of the zero of B). The 
maximum of B is outside the range of the plot, i?max > 500. 



Both types of solutions have the property that the Riemann tensor and Ricci 
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Figure 24: Singular shell case. Function A(r) for A = 1 and k = 0.1. Shooting 
from the compacton boundary for the initial values nis = 0.5 and R = 2. The 
resulting singular boundary is at tq = 0.9928 (position of the zero of A). 
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Figure 25: Singular shell case. Product A{r)B{r) for A = 1 and k = 0.1. Shoot- 
ing from the compacton boundary for the initial values = 0.5 and R = 2. The 
resulting singular boundary is at tq = 0.9928 (position of the zero of AB). 



tensor are singular at r = 0. Singularities show their presence in a most physically 
relevant way in invariants which may be derived from the Riemann and Ricci 
tensors. A further simplification occurs for invariants derivable form the Ricci 
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tensor alone, because then the Einstein equations may be used. Concretely we 
find for the Ricci scalar 

n = -K^Ti; = Va02 (58) 

so the Ricci scalar is, in fact, regular at the origin or the singular boundary, because 
the scalar field (p is. The Ricci tensor squared may be expressed like 

7^^,7^'^- = «4 ^j.^^j.,u ^ ^r^^Y) = + (59) 

where we used 



T,.T^''=-^+iXY- (60) 

Here, the first term is singular at r = 0, because 0' is regular and non-zero there 
whereas B goes to zero like r (large compacton), (small compacton), or (r — 
To) (singular shell), respectively. Therefore, the Ricci tensor squared is already 
singular. More complicated invariants like, e.g., the Kretschmann invariant, are 
singular, as well. 

One possible question to be asked in connection with these singularities is 
whether they really belong to the space time manifold, that is, whether a freely 
falling particle may reach them in finite proper time. We will find in all three 
cases that this is the case, i.e. there exist geodesies which hit the singularities in 
finite proper time. For this purpose we have to study the geodesic equation 

^''(^) + = 0. (61) 

Here is the Christoffel connection and r is the proper time. We shall restrict 
to radial geodesies, which is sufficient for our purpose. In this case the geodesic 
equations reduce to 

r + r;/ + 2r;,ir + r;/' = o (62) 

where the nonzero components of the Christoffel connection are 

1 1 A' 

1 1 A' 

r„ = ^g"drgrr = l^ (63) 
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We are interested in the geodesic motion near the singularity, therefore it is enough 
for our purposes to restrict to the leading behaviour of A and B. In the large 
compacton case , this leading behaviour is ^4 ~ (a_i/r) and B ~ bir and we get 
the geodesic equations 



i --if = 
r 



where 7 = (a_i/6i). The equation for t has the solution 

t^ci r dT'r{T') + C2 (65) 



and the equation for r becomes 



C3 , 1 .2 



r - — + —r = (66) 
r 2r 

where C3 = (70^/2) > 0. Under the transformation r — the last equation turns 
into 

/ = (67) 

This last equation has an easy interpretation in terms of an equivalent mechanical 
system. It is just the e.o.m. of a nonrelativistic particle in one dimension in the 
repulsive external potential U{f) ~ — . The center / = of the potential may 
be reached despite its repulsive nature because of the positive power (2/3) in the 
distance law. However, not all radial geodesies hit the center r — 0, because in 
the equivalent mechanical problem the particle needs a sufficient initial velocity. 

The small compacton case may be analysed in an equivalent manner. The 
leading behaviour is >1 ~ a2r^ and B ~ 62r^ and we get the geodesic equations 

.. 2 . 

t+-tf = 
r 

f+^P + ^f'^ = (68) 

where 7 = (02/62)- The equation for t has the solution 

r dr' 

t = ci -^r^ + C2 (69) 



To 



r2(T') 
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and the equation for r becomes 



f--^- (VI) 



r + ^ + -r^ = (70) 

where C3 = (70^/2) > 0. Under the transformation r = the last equation turns 
into 

2C3 

p 

This is the e.o.m. of a nonrelativistic one-dimensional particle in an attractive 1// 
potential. The particle will hit the center / = unless it has a sufficiently large 
(outward directed) initial escape velocity. 

Finally, for the singular shell it is useful to introduce the variable u — r — vq. 
Then the leading behaviour is A ~ a^u^ and B ~ biu and we get the geodesic 
equations 

.. 3 . 

t + -f/u = 
u 

+ — +7r^ = (72) 
2 2u 

where 7 = (03/61). The equation for t has the solution 

r dr' 

t^ci —— + C2 (73) 



and the equation for u becomes 



C3 , 1 .2 



+ — K^ = (74) 

where C3 = (870^/2) > 0. Under the transformation u = the last equation 
turns into 

This is the e.o.m. of a nonrelativistic one-dimensional particle in an attractive 
1/P potential. The particle will hit the center / = unless it has a sufficiently 
large (outward directed) initial escape velocity. 

We conclude that in all cases there exist radial geodesies which hit the singu- 
larities at r = or r = ro, respectively, therefore these singularities belong to the 
corresponding space time manifolds. 
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The singular shells, however, differ from the compact balls in one essential 
aspect, in that their singular inner boundary is, at the same time, the locus of a 
Killing horizon and, therefore, invisible for an outside observer. In a first step, let 
us consider the equation for a light-like radial geodesic. We find 

iB 

ds^ = = -Adt^ + Bdr^ ^ dt = \ —dr (76) 

\ A 

and with B/Ar^ (c/ (r — ro))^ near r = ro we find the solution 

r(i) - ro = exp (77) 

where tg is an integration constant. Obviously, r = tq requires t = — oo, so a 
light ray cannot escape from the singular surface. A slightly more rigorous and 
less coordinate dependent derivation uses the concept of a Killing horizon. Let us 
assume the existence of a Killing vector ^ and consider the set of points (hyper- 
surface) Hq where is null, i.e., = ('CiO = 0. A Killing horizon if is a 
connected component H E Hq which is, at the same time, a null hyper- surface, 
i.e., M = [dN, dN) = 0. So let us demonstrate that the inner boundary of the 
singular shell is indeed a Killing horizon. For the Killing vector we choose the 
generator of time translations ^ = (d/dt) = dt, as is obvious for a static space- 
time. For N we find 

N = (e, = 9^.udx'' ® dx'idt, dt) = -A. (78) 

As near r = tq A behaves like A ~ a3(r — vq)^, it holds that N{ro) — 0. For M 
we get with dN = —A'dr 

M = {dN, dN) = A'^g^'^d^ ® d^{dr, dr) = —. (79) 

B 

With the leading behaviour of A near ro like above and B ~ 6i(r — ro), we 
easily conclude that indeed M(ro) = 0, so the singular boundary is a Killing 
horizon. We remind the reader that for static space-times every event horizon is 
a Killing horizon, which makes our result all the more interesting. The singular 
shell solutions we found are, in fact, quite similar to ordinary black holes in this 
respect. 
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4 Discussion 



It has been the main purpose of the present article to investigate in detail the prop- 
erties of compact boson stars (compact balls or shells minimally coupled to grav- 
ity) for a theory with a non-standard kinetic term in three plus one dimensions. 
We chose the simplest possible theory for this purpose, with a purely quartic ki- 
netic term and the simple quadratic potential, because we wanted to pursue the 
analytical investigation as far as possible in order to understand also more generic 
and qualitative features, in addition to explicit numerical calculations. This theory 
has a one-parameter family of compact ball solutions already in the case without 
gravity. The compact balls in the theory without gravity turn out to be solutions 
of the weak type, which solve the field equations everywhere except at the ori- 
gin. We remark that this fact is not a problem in the case at hand, because the 
set of weak solutions is the appropriate solution set for variational problems of 
this type, in any case. Alternatively, the compact ball solutions may be extended 
to solutions in the whole space by introducing a delta function source term at the 
origin. It turned out that these solutions may have arbitrary size and energy. 

In the case with gravity, we found two different types of solutions, both of 
which are of the strong type, i.e., they solve the field equations everywhere. For 
compact balls, we found that the scalar field itself is behaving well at the origin 
(it has a power series expansion there). The metric coefficients do have a power 
series expansion at the origin, as well, (with the exception of a r^^ leading term in 
one specific case), but nevertheless give rise to a singularity at r = 0. Specifically, 
some higher invariants of the Riemann tensor, like the Ricci tensor squared, or 
the Kretschmann invariant, become singular at this point. This singularity is not 
shielded by any horizon, so it is a naked singularity. In spite of the naked sin- 
gularity at the center, asymptotically these solutions behave like Schwarzschild 
solutions, so they have finite Schwarzschild mass (or asymptotic, ADM mass) 
nis. In one class of solutions (the large compactons), this ADM mass may take 
on negative values, in spite of the positive definite energy density. The theorems 
excluding this type of behaviour are evaded by the presence of naked singularities 
in the case at hand. 

In addition, we found another type of solutions in the theory with gravity, 
namely the singular shells where a singularity already forms at a finite, nonzero 
value of the radial coordinate. This singularity is a surface in space with finite 
area, as may be checked easily by inspecting the corresponding metric (the angu- 
lar part of the metric is not multiplied by any suppressing factor at the position 
r = ro of the singularity). Further, this singularity is a genuine singular boundary 
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of space-time, because space-time cannot be continued beyond this boundary. At 
the same time, this singular boundary is a Killing horizon. In addition, the sin- 
gular shells asymptotically for large radius are Schwarzschild, so they resemble 
ordinary Schwarzschild black holes in many respects. Certainly they are black 
(i.e., no radiation may escape from the horizon), although the curvature near the 
horizon is different from the Schwarzschild black hole case. It is probably more 
a question of parlance whether these objects should be called genuine black holes 
- thereby providing a counter-example to the scalar no hair conjecture - or inter- 
preted as different object which, nevertheless, share many properties with ordinary 
black holes. In any case, the existence of these singular shell solutions is interest- 
ing. On the other hand, ordinary Schwarzschild black hole solutions (where the 
curvature near the horizon behaves like in the Schwarzschild case) do not exist in 
our model, except for the trivial case where the metric is exactly Schwarzschild 
everywhere and the scalar field takes its vacuum value in all space-time. 

At this point, an obvious question to ask is whether theories of the type pre- 
sented here may be of some relevance in astrophysical or cosmological contexts. 
The theory presented here is the simplest possible case of the class of theories 
having a non-standard kinetic term and allowing for compact boson stars. There- 
fore, an attempt for direct applications would most likely be premature, and we 
shall focus, instead, on more generic issues. Specifically, we want to comment 
on some features of these theories which might make them interesting for as- 
trophysical considerations. A first feature we already mentioned is the fact that 
non-overlapping multi compact ball configurations do not interact at all in the 
case without gravity, so the self-gravitating ones only interact via the universal 
gravitational interaction. A further feature of compact solitons in K field theories 
in general is the fact that linear fluctuations are possible only inside the compact 
solitons, whereas they are completely suppressed in the vacuum. This fact is easy 
to understand. Unlike for the standard kinetic term, in the case of a non-standard 
K theory the wave operator acting on the fluctuation field 50 is multiplied by some 
power of X = d'^(t)d^(i), where X has to be evaluated at the background field (i.e., 
the compacton). But X is identically zero outside the compacton (that is, in the 
region where the background field takes its vacuum value), which is the reason for 
the complete suppression of linear fluctuations in that region. This implies that a 
compacton in K field theories cannot emit particle type radiation, which makes it a 
rather stable object. These two facts - that compact boson stars only interact grav- 
itationally, and that particle-like excitations in vacuum do not exist - might make 
compact boson stars of the type discussed in this paper interesting candidates for 
dark matter. The second fact may, further, give rise to the following speculations. 
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Let us assume for a moment that a universe with a scalar K field theory was a 
reasonably good effective description in a certain (early) stage of the evolution of 
the universe. Then linear fluctuations of the K field which may act as seeds for 
matter formation were naturally present in the compacton regions, but absent in 
the vacuum regions. The late time remnants of the former would be regions of the 
universe with a higher matter concentration, whereas the latter might have evolved 
into voids (low matter density regions). 

Before more convincing conclusions can be drawn with respect to possible 
applications of compact boson stars, one key issue to understand is which features 
of the theory presented in this article are generic and which are just artifacts of the 
specific model. One important feature in this respect obviously is the presence of 
naked singularities in the compact ball type solutions of the model studied in this 
paper. Such naked singularities of the point-like type are already present in finite 
mass solutions (boson stars) for standard self-gravitating scalar fields f29\ - [|32ll . 
Further, it is known that the gravitational collapse of a standard scalar field may 
terminate in a naked singularity (391, iBDl- Nevertheless, the presence of naked 
singularities is frequently viewed as an unwanted or even unphysical feature and, 
in any case, the generic or non-generic nature of naked singularities in compact 
balls is an interesting question. We cannot give a completely mathematically rig- 
orous answer to this question, but shall try to partially answer it relying on our 
experience and some recent results in the literature. On the one hand, for a static, 
radially symmetric scalar field, it seems to be a typical feature that non-singular 
metric functions A, B at the origin r = require that the scalar field takes its 
vacuum value there. Further, these local solutions typically are isolated solutions 
(i.e., without free integration constants) which makes it difficult to connect them 
to finite mass (i.e., Schwarzschild) asymptotic solutions. In our specific model, 
we could exclude this possibility, whereas in more general theories it might still 
be possible to connect a regular center with an asymptotic Schwarzschild solution 
provided that i) the scalar field has more than one vacuum and ii) the scalar field 
Lagrangian contains sufficiently many coupling constants which may be varied in 
order to connect the two local solutions near zero and near infinity. 

There exists, however, a class of slightly different but related theories where 
the appearance of compact balls with a regular center seems to be a more generic 
feature, namely the class of self-gravitating compact Q balls. A Q ball is a com- 
plex scalar field ^ restricted to the ansatz tjj = exp{iujt)(f)(x), such that the re- 
sulting static equations for the real function (p are quite similar to the static field 
equations of a real scalar field, both with and without gravitation. Q balls (without 
gravitation) were introduced by S. Coleman in [|41l and have been widely inves- 
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tigated since then. Systems of self-gravitating Q balls are frequently subsumed 
under the general notion of boson stars. For more details we refer to the review 
articles [36J - fi38il . Compact Q balls (without gravity) for a complex scalar the- 
ory with a standard kinetic term but with a V-shaped potential have recently been 
found in [|6l and in Q (the second paper included also an electromagnetic inter- 
action of the complex scalar, and allows both for ball-shaped and shell-shaped 
solutions). The gravitational counterpart of the model studied in [71 was investi- 
gated in [28J, mainly using numerical methods, and the existence of non-singular 
solutions (both of the ball type and of the shell type) was established. Therefore, 
the existence of regular self-gravitating compact solutions seems to be a more 
generic feature for Q balls, analogously to the case of standard (non-compact) Q 
balls. 

In any case, we think that the further study of the interaction between compact 
objects (like compact solitons, Q balls, etc.) and gravitation is an interesting en- 
deavour which may lead to new insights and unexpected phenomena, as well as 
some applications to astrophysics and cosmology. One future line of investigation 
will certainly be the study of Q balls. Another promising direction is the analysis 
of the symmetries and conservation laws of the corresponding theories, e.g., along 
the lines of [i42l - P31 . which may lead to a better understanding of more general 
solutions as well as their stability. These issues are under current investigation. 
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